Kovacs effect in facilitated spin models of strong and fragile glasses 
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We investigate the Kovacs (or crossover) effect in faciiitated /-spin modeis of gfassy dynamics. 
Although the Kovacs hump shows a behavior qualitatively similar for all cases we have examined 
(irrespective of the facilitation parameter / and the spatial dimension d) , we find that the dependence 
of the Kovacs peak time on the temperature of the second quench allows to distinguish among 
different microscopic mechanisms responsible for the glassy relaxation (e.g. cooperative vs defect 
diffusion). We also analyze the inherent structure dynamics underlying the Kovacs protocol, and 
find that the class of facilitated spin models with d > 1 and / > 1 shows features resembling those 
obtained recently in a realistic model of fragile glass forming liquid. 
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Introduction 

There is much interest in understanding how the out 
of equilibrium properties of glassy materials can be de- 
scribed within a statistical mechanics framework. At low 
temperatures, the exceedingly long relaxation time pre- 
vents such materials from reaching thermal equilibrium 
on the accessible time scales, and makes hard to establish 
a connection with the thermodynamic phase behavior . 
In particular, several questions about the nature of the 
glass state are still unanswered: for example, does it rep- 
resent a genuine thermodynamic phase of the matter or 
a purely kinetic phenomenon? How many "Active" ther- 
modynamic parameters are needed for characterizing the 
amorphous state? These and other questions are relevant 
not only for glasses but for every system with slow dy- 
namics 2], such as dense granular materials and colloids, 
disordered dielectrics and electron glasses, etc. 

Fifty years ago, Kovacs designed an experiment in 
which the low temperature behavior of glasses is probed 
by means of the following protocol (see also 0, 0)- 
A sample is brought from high to low temperature Xi, 
and left to age for a time t w until its volume v attains 
the value it would have in thermal equilibrium at a tem- 
perature T% > Ti, v(t w ) — w cq (T2). The sample is then 
suddenly heated to T2 and the ensuing volume relaxation 
is observed. Even though the average volume is already 
at its equilibrium value, it turns out that in order to fully 
equilibrate, the system first increases its volume, then 
passes through a maximum that depends on the ther- 
mal history of the system, and eventually recovers the 
equilibrium. The non-monotonic behavior generally re- 
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suits from the superposition of two antagonist processes: 
the fast modes - equilibrated at T%, tend to increase the 
volume, while the slow modes - which are still relaxing to- 
wards T\ , tend to decrease it Q . The sample is therefore 
formed by regions having an average local density larger 
or smaller than the equilibrium one, which respond on 
different time scales to a change in temperature. The 
experiment thus provides an indirect evidence that the 
low temperature dynamics of glasses is spatially hetero- 
geneous nnn. 

The Kovacs effect was first observed in polymeric 
glasses and in statistical mechanical systems as sim- 
ple as the ferromagnetic Ising chain 0. There is an in- 
creasing interest in this phenomenon in connection with 
other kind of memory effects investigated in disordered 
magnets, and recent statistical mechanical approaches to 
glassy and granular matter. In these contexts, it has 
been observed that the Kovacs effect occurs in systems 
with quasi-long range order pfl, spin glasses |g, realis- 
tic model of molecular liquid (ortho-terphenyl) traps 
models [l(| , kinetically constrained spin chains , and 
in simple models of granular compaction, such as the 
parking lot model ^3 ■ More recently, the effect has also 
been studied in the spherical p-spin glass |l3j . 

In this paper we study the Kovacs effect in facilitated 
/-spin models of glassy dynamics in spatial dimension d. 
Our interest here is mainly dictated by the possibility of 
using the Kovacs effect to identify the microscopic mech- 
anism responsible for the slow relaxation (e.g., strong vs. 
fragile glass behavior 0]), or ' more ambitiously, to dis- 
criminate between situations in which an arrest transi- 
tion occurs from others in which it does not. We also 
investigate the inherent structure dynamics during the 
Kovacs protocol and show that models in spatial dimen- 
sion d > 1 and facilitation parameter / > 1 are able to 
account for a subtle effect recently observed in molecular 
dynamics simulation of a realistic model of fragile glass 
forming liquid 



2 



Facilitated spin models 

The model we consider is model first introduced by 
Fredrickson and Andersen |l5j . with the purpose of de- 
scribing the cooperative relaxation in highly viscous liq- 
uids. The spins are non interacting and the Hamiltonian 
density is simply 
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(1) 



where the s, = 0, 1 are spin variables, and the index i 
runs over the sites of a lattice with coordination number 
c. One possible interpretation of the spin variable is that 
its state represents, after a suitable coarse-graining over 
molecular length and time scales, mesoscopic regions of 
the liquid with more or less density. What makes the 
model non trivial is the dynamics, which is defined by 
the following facilitation rule: Sj can flip with a non-zero 
rate if and only if at least / of its c neighboring spins are 
up, that is 



hi = Y S k > f ■ 



(2) 



The transition rates satisfy detailed balance and are pro- 
portional to the Metropolis factor 



W(s % -v 1 - Si ) = (hi - f) min 1, e 



(3) 



where /3 is the inverse temperature (3 = l/k B T, and the 
Heaviside function © encodes the facilitation mechanism. 
It turns out that at low enough temperature the sys- 
tem dynamics is unable to relax to the equilibrium dis- 
tribution within the simulation time window. In fact, 
as the density of up spins becomes smaller and smaller, 
the facilitation rule is hardly satisfied and consequently 
the spin relaxation becomes slower and slower. For this 
reason the model exhibits several glassy features, such 
as non-exponential relaxation, history dependent prop- 
erties and hetero gene ous distribution of local relaxation 
times [la lla ll3 ■ F° r a review see Refs. |l9j. This 
type of facilitated spin model is reminiscent of the ki- 
netically constrained lattice-gas introduced by Kob and 
Andersen 20] . However, the facilitated dynamics is not 
conserved and the constraint is weaker (as it involves 
only the local field hi of the spin which attempts to flip). 
Numerical investigations of two and three dimensional 
systems (with different values of / > 1) suggested that 
the characteristic equilibration time at low temperature 
diverges as a Vogel-Fulcher law [lTj, although early ar- 
guments shown that the s yste m is ergodic at any tem- 
perature larger than zero ]WL l2l| . More generally, the 
equilibration time has a non-Arrhenius behavior typical 
of fragile glasses for d > 1 and / > 1. The case / = 1 
instead corresponds to the defect diffusion dynamics and 
leads to the Arrhenius behavior, logi cq ~ 1/T, charac- 
terizing strong glasses. However, Id models, such as the 



East model |23j, |2J| which have directed (or asymmetric) 
constraints, may also exhibit a non-Arrhenius behavior. 

One reason of interest in the facilitated spin models 
is that they present many glassy features in spite of the 
absence of a true glass transition in finite dimension pH 
I22T ]. while on the Bethe lattice they appear to have an 
ergodic/non-ergodic transitions of the kind predicted by 
mode-coupling theory |25|, l2(il l27j . In the following we 
shall be interested in reproducing the Kovacs experiment 
in facilitated spin models defined on finite dimensional 
lattices with the purpose of comparing the behavior for 
different value of the facilitation parameter. 



Kovacs effect 

To numerically reproduce the Kovacs experiment 
an equilibrium configuration at high temperature is 
quenched at low temperature T% and left to age for a 
waiting time t w , until the energy density, e, attains the 
value it would have in thermal equilibrium at tempera- 
ture T 2 > Ti, e(t w ) = e oq (T 2 ). According to Eq. JTJ the 
equilibrium energy density reads: 
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The energy density therefore is the natural observable 
that plays in our model the role of the volume in the Ko- 
vacs protocol. This quenching procedure is repeated for 
several temperatures T\. Each configuration obtained at 
time t w such that e(i w ) = e oq (T2), represents a distinct 
glassy state with different thermal history but identical 
energy. Such configurations are then quickly heated at 
temperature T 2 > Ti, and the reduced energy Av, de- 
fined as 



A«(Mw) 



e(t) - e(t w ) 
e(t w ) 



(5) 



is recorded as a function of time t (t > t w ). We use a 
fast algorithm that keeps track of the list of free (un- 
constrained) spins, from which the attempted spin-flips 
are chosen with the usual Metropolis algorithm. Accord- 
ingly, the Monte Carlo time sweep is rescaled with the 
number of mobile spins. The largest number of spins we 
consider is N — 2 21 . The lattice systems have periodic 
boundary conditions. 

In Fig. ^ (top) we show the reduced energy relaxation 
during a numerical experiment for a 3d facilitated spin 
system with / = 3, final temperature T2 = 1.0 and 
several initial quenching temperatures T\. Qualitatively 
similar results are obtained for d — 3 and smaller facili- 
tation parameters (/ = 1, 2), for the square lattice with 
/ = 1, 2, as well as on a Bethe lattice (a random graph 
with fixed connectivity) where the glass transition occurs 
at a finite temperature for / > 1 27]. The plot exhibits 
the typical Kovacs humpvery similar to those obtained in 
many other systems |SB II H ES IH EEl , the hump 
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FIG. 1: Kovacs effect in a facilitated spin model on a cubic lat- 
tice with facilitation parameter / = 3 (system size V = 128 3 ). 
The system relaxes at temperature Ti until it attains the en- 
ergy e(i w ) = e cq (T2); afterwards the temperature is changed 
from Ti to T 2 = 1.0 (top) or 0.8 (bottom). The full line repre- 
sents the energy relaxation after a quench at T2. Qualitatively 
similar results are found for different dimensions, facilitations 
and geometries. 



height increasing with AT = T2 — T\ and shifting to 
smaller times as T\ decreases. By decreasing AT one ob- 
serves (see Fig.^ bottom, where T 2 = 0.8) the clear two- 
step behavior of the energy relaxation, as recently found 
in the Id case The fast, non activated processes 

due to the excitations of spins neighboring facilitated 
ones, contribute to the formation of a plateau starting 
on timescales t — i w ~ 0{1). On longer timescales, 
due to activated processes, the energy increases up to a 
maximum at time t — t w = t K , decreasing afterwards un- 
til the equilibrium is eventually attained. When rescaling 
time by a suitable characteristic time, t/t*, and Av by 
the maximum height, Aw max , a good collapse of the late 
time regime is obtained, while deviations occur for the 
non activated part of the hump, as can be seen in Fig. EI 
Usually t* has to be chosen much larger than the time 
of the maximum, i K , so that the contribution of the fast, 
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FIG. 2: Scaling of the data of Fig. □ for T 2 = 0.8. Both 
time and height are rescaled, respectively by t* (see text) and 
Au max . Data collapse is good except for the initial part where 
scaling is not expected to hold @- 
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FIG. 3: Kovacs effect in a facilitated spin model on a cubic 
lattice with / = 3 (system size V = 128 3 ). The protocol here 
is different from fig. the system relaxes at temperature T\ 
for a time t w = 1000 and then the temperature is changed 
from Ti to T 2 = 1. See also Ref. II. 



non scaling modes vanishes Here, instead, as a con- 
sequence of the clear separation of time scales, t* is very 
close to t K . 

Another possible way to investigate the Kovacs effect 
is to fix the time spent at temperature Ti and then to 
change the temperature to Ti. An example is shown in 
Fig. [31 for d = f = 3. This protocol has been used in 
Ref. [9( . One can notice that the results are qualitatively 
similar to those obtained with the previous protocol. 

Finally, we must remark that, in order to remove the 
trivial part of the Kovacs effect, one should consider 
the initial time slightly after the jump as discussed in 
Ref. jlfj. However, for the class of models considered 
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here the contribution of fast modes, lasting for 0(1) MCs, 
is clearly resolved due the presence of an intermediate 
plateau. Moreover, this initial increase can be removed 
by looking at the inherent structures dynamics. We also 
checked that within our protocol there is no Kovacs effect 
in the absence of kinetic constraints, confirming that the 
observed effects are not an artifact of our procedure. 



Temperature dependence of the Kovacs peak 

The results we have have obtained for different values 
of the facilitation parameter /, would lead to the con- 
clusion that the Kovacs effect is unable to discriminate 
among the several microscopic mechanisms of slow re- 
laxation or, even worse, between situations in which a 
dynamical arrest transition occurs from others in which 
it does not. As we show below, this is not the case: a 
more careful analyzes shows clear differences between all 
these cases. In particular, we have analyzed the behav- 
ior of the Kovacs peak time t K as a function of the final 
temperature T2 at fixed Ti, making a comparison with 
recent estimations of the equilibration times coming from 
independent methods. 

For d = 1 there are simple arguments 1 1 1 1 | that cor- 
rectly describe the dependence of both t w and t K on T\ 
and T 2 . The average distance between up spins is pro- 
portional to the inverse of their concentration, e _1 (T) ~ 
exp(l/T), for low T, while the effective diffusion rate for 
these spins is T(T) ~ exp(— 1/T). Thus, after the initial 
quench to Ti, the time t w to attain the equilibrium con- 
centration (at T2) of up spins, e eq (T 2 ) ~ exp(— I/T2), 
is t w (Ti,T 2 ) ~ r-^TiJe-^Ta) ~ exp(l/T 1 + 2/T 2 ). 
When T\ = Ti we obtain the equilibration time af- 
ter a quench, t oq ~ exp(3/T2). A similar argument, 
for the East model, leads to i oq ~ exp(l/T 2 2 log 2) and 
i w ~ exp(l/T]T 2 log2). In both cases, it turns out that, 
as a function of T 2 , t K ~ i oq . However, it is not clear 
whether simple arguments exist for other dimensions and 
facilitation parameters, although the correspondence be- 
tween i K and i oq seems to remain valid, as is shown below. 

Fig. 0] (top) shows t K vs X2 for the defect diffusion 
dynamics, / = 1, and spatial dimension d < 4. In 
all cases, t K shows an Arrhcnius behavior in the limit 
T2 — > T x + and scales with T2 as the equilibration re- 
laxation time of the system. For d = 1 we recover the 
result of Ref. [TJ, with t K ~ exp(3/T 2 ). For higher di- 
mensions, a generalization of the previous argument [l9| 
gives t oq ~ exp[(l + 2/<f)/T 2 ], and seems not to be correct 
(dashed lines in figure). Indeed, for d = 2, the result t K ~ 
exp(2.3/T2) is compatible with many, independent nu- 
merical estimates of t eq , for example, from the integrated 
correlation function and the dynamical susceptibility psj 
and from the distribution of persistence times [^j . This 
value is in the same universality class of directed perco- 
lation (DP) I3C| . For d = 3 and 4, a renormalization 
group analysis [3l| gives i cq ~ exp[(2 + e/12)/T 2 ], with 
e = 4 — d, again in the same universality class of DP. 
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FIG. 4: Behavior of t function of T a -1 . Top: / = 1 

and Ti = 0.2. In all cases the behavior, for T2 — > 
is Arrhenius. For comparison we show several predictions 
available in the literature. The solid lines correspond to: 
t K ~ exp(2.3/T 2 ) HI |H i3 (d = 2) and t K ~ exp[(2 + 
e/12)/T 2 ] H3) with e = 4 - d (d = 3 and 4). The dashed 
lines correspond to the mean field result [ill . Il9| : t K ~ 
exp[(l + 2/d)/T2\. Bottom: d = f > 1, showing super Arrhc- 
nius behavior (see text). Notice the iterated exponential in 
the abscissa: for example, exp' 2 '(a;) = exp(exp(a;)). 



In Fig. ^] (bottom) we show the results for the coop- 
erative diffusion dynamics with d = f > 1. In this case 
the systems exhibit a fragile glass behavior. Accordingly, 
we checked the Kovacs peak against an iterated expo- 
nential function t K ~ exp(exp(6 2 /T)) for d = f = 2, and 
t K - exp(exp(exp(6 3 /T))) ford = / = 3]lllllEal23- 
Again, our estimations agree pretty well, within the nu- 
merical accuracy, with those of the equilibrium relaxation 
time, obtained independently from the integral of the cor- 
relation function and the peak of the dynamical suscep- 
tibility [23. For the East model, whose dynamics is also 
cooperative, again there is a correspondence with the re- 
laxation time of the system, A first conclusion that we 
can draw from this analysis is that the Kovacs protocol is 
rather sensitive to the strong/fragile relaxation behavior 
and seems to provide an alternative, independent method 
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FIG. 5: Energy density of inherent structures during the Ko- 
vacs protocol for d = f = 3 and several values of Ti. The 
lattice side is L = 64 and T2 = 0.8. The arrows show two 
ISs (with energy ei S = 0.2179 and 0.2158) obtained by a zero- 
temperature descent from configurations having identical en- 
ergy e = 0.24 before and after the peak. 



even if the energies at time t w are the same for all Ti, 
e(t w ) — e cq (T2), the corresponding inherent structures 
arc different: the lower is Ti, the higher is e IS (t w ). For 
times up to 0(1), the IS evolution is confined in the same 
"energy basin" since the fast thermal excitations are re- 
moved by the zero-temperature descent. On longer time 
scales, the slow, activated processes start to play a role, 
the system is able to overcome the energy barriers and 
leave the initial basin, given origin to the hump. After 
the hump, the curves merge with the one obtained from 
a direct quench to T2 (thick line in Fig.|SJ) and eventually 
reach a value different (lower) than the initial one, that 
only depends on the equilibrium energy, e IS (e eq ) (e.g., for 
/ = 1, e IS (e oq ) = e e ^l — e eq ) d [36|). Thus, in analogy with 
the results of ref. here too the IS explored far from 
equilibrium differ in depth from those sampled close to 
equilibrium. This difference increases with T2 — Ti and is 
due to a different number of excitations, or equivalently, 
to a different realization of spatial heterogeneities, and 
not to topological properties of distinct regions of the 
phase space. 



to obtain the equilibration time of the system. 

Some extra information can also be obtained from t w , 
the time spent at T±. In all cases we have considered here, 
the hump time, t K , and t w {Ti,T 2 ), are related through a 
power-law t K ~ i^, where b is temperature independent 
if d = /, and symmetric facilitation. On the other hand, 
for d > f = 1 and in the East model, b ~ I/T2 and 
b ~ T2, respectively. In domain growth models one finds: 
< K - <w /2 (2d XY 7]) and t K ~ t z J (1+z) for the Ising 
model |Hj. For d = f = 1 we find b ~ 0.68, close to 
the expected 2/3 because in Id the dynamics proceeds as 
a domain growth process with dynamic exponent z = 2 
(we also obtain 2/3 from the expressions for i K and t w 
derived in Ref. 11]). On the other hand, for d = f = 2 
and 3 we get b ~ 0.89 and 1.05, respectively, close to 
unity, as in the p-spin-glass model ^3j, where b ~ 0.9. 



Inherent Structure Dynamics 

Finally, in order to make contact with recent results ob- 
tained by molecular dynamics simulation, Ref. B, w e in- 
vestigated the inherent structure (IS) dynamics |34( dur- 
ing the Kovacs protocol. For spin models with stochastic 
dynamics such as those considered here, an IS is a blocked 
configuration, i.e. a configuration in which only isolated 
defects are present and no spin-flip is possible. The IS dy- 
namics is obtained as follows: at every time t a copy of 
the system is made and evolved with a zero-temperature 
sequential dynamics until it reaches a configuration in 
which no further spin-flip is possible. The energy e IS 
of this IS is recorded and the dynamics of the original 
systems is resumed, this procedure being repeated along 
the system evolution. Fig. [S] shows a typical example of 
IS dynamics for d = f = 3. First, one observes that 



Once the system leaves the initial IS, two different be- 
haviors are possible when one considers a couple of ISs 
descending from configurations having the same energy 
before and after the peak. While in the simplest case the 
two ISs may have the same energy there could be a more 
complicated case for which this is not true. The latter 
case has been observed recently by Mossa and Sciortino 
in molecular dynamics simulation of orto-terphenyl |9(, 
a typical fragile glass forming liquid. We find that for 
/ > 1 there is an IS energy-shift similar to that observed 
in Ref. (see Fig. top, and the arrows in Fig. [SJ . 
For 1-spin facilitated dynamics instead, the ISs have the 
same energy (see Fig. [51 bottom, for an example) . The 
reason of this distinct behavior does not seem directly 
related to the form of the relaxation time, but rather to 
the value of the facilitation parameter /. In fact, we find 
that in the East model, whose relaxation time is super- 
Arrhenius, no IS energy-shift is observed either 37] . Our 
results seem to suggest that only facilitated spin models 
with d > 1 and / > 1 are able to reproduce the subtle 
effect observed in Ref. 



In Fig. we also show the IS sampled in equilibrium 
(thick lines) and, interestingly, in all cases, during the 
whole protocol the system explores regions never vis- 
ited while in equilibrium: only asymptotically the points 
merges with the equilibrium curve. However, for smaller 
temperatures, the separation between the lower branch 
and the equilibrium curve diminishes, as well as the dif- 
ference between the branches. Notice that in Ref. [jj, 
soon after the hump, the system starts to sample config- 
urations from a region of the phase space that is explored 
in equilibrium (but at different temperatures). 
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FIG. 6: Energy density of inherent structures during the Ko- 
vacs protocol versus the energy density. Top: d = / = 3 
for Ti = 0.6 and T2 = 0.8. The thick line is the equilibrium 
value, and both curves merge only at the end of the Kovacs 
protocol. Bottom: d = / = 1 for Ti = 0.4 and T 2 = 0.7. 
Notice that in this case, there is no difference between config- 
urations with the same energy on both sides of the hump. For 
late times the curve merges with the equilibrium one (thick 
line), eg = e cq (l - e cq ) d . 



Conclusions 

In summary, we have studied the Kovacs effect in a 
class of spin models with facilitated dynamics which ex- 



hibits slow relaxation at low temperature. Our results 
show that the Kovacs effect is qualitatively similar in 
all cases examined, irrespective of the spatial dimension 
the facilitation parameter, and the presence or not of 
a dynamic transition. However, in the class of facili- 
tated models considered here one may distinguish be- 
tween these situations by looking at temperature (of the 
final quench, T2) dependence of the Kovacs peak time, 
which scales in the same way as the equilibrium relax- 
ation time of the system at that temperature. Thus, the 
Kovacs protocol provides an independent way to gain in- 
formation about the nature of the relaxation dynamics, 
at least for this class of models. It would be interest- 
ing to understand to what extent such a correspondence 
holds in more realistic systems, or in systems where the 
structural relaxation is related with the topography of 
the underlying potential energy landscape. 

A related question concerns the inherent structures vis- 
ited while the system follows the Kovacs protocol. Gen- 
erally, they show the non trivial, activated character of 
the Kovacs effect, as well as the role played by spatial 
heterogeneities. In contrast with Ref. the IS vis- 
ited by the facilitated models considered here are not the 
equilibrium ones, which are only attained asymptotically. 
Moreover, the East and facilitated 1-spin models are not 
able to reproduce the subtle IS energy-shift observed in 
Ref. @, while such an effect occurs in facilitated spin 
models with d > 1 and / > 1. This suggests that the 
inherent structure dynamics underlying the Kovacs pro- 
tocol and the IS energy-shift might be used as a tool for 
discriminating between distinct microscopic mechanisms 
of dynamic facilitation. It would be also interesting to 
check by molecular dynamics simulation whether simi- 
lar results hold for realistic models of strong and fragile 
glasses. This could provide a more stringent test of dif- 
ferent scenarios of glassy dynamics. 
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